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I. INTRODUCTORY REMARKS 

In the 60-es and 70-es, when the basic concepts 
of dynamical symmetries have been formulated and 
elaborated^— only few physical realizations of these sym- 
metries could be found in the realm of existing quantum 
mechanical objects (see 7 -^ for a review). On may men- 
tion the hydrogen atom^r— and the harmonic oscillator in 
various spatial dimensions^— as the systems for which 
the study of their dynamical symmetries revealed addi- 
tional facets of excitation spectra and response to ex- 
ternal fields. Rapid progress in the nanotechnology and 
nanophysics during two recent decades significantly ex- 
tended the field of applicability of these concepts and 
enriched the theory with some new ideas. 

Contemporary nanophysics^— deals with the artifi- 
cial structures which consist of finite number of electrons 
confined within a tiny region of space, where the energy 
spectrum of electrons is discrete. As a result, such ob- 
jects can be treated as " zero-dimensional" artificial atoms 
or molecules with spatially quantized discrete states, well 
defined symmetry and controllable electron occupation. 
Besides, modern technologies allow fabrication of devices 
where a "natural" atom or molecule is spatially isolated 
from the rest part of a device, so that the physical prop- 
erties of an individual atom or atomic cluster may be 
studied experimentally. 

In this paper we analyze the dynamical symmetries 
which arise when a group theoretical approach is used 
in the description of a contact between a few electron 
nanosystem S with definite symmetry Gg and a macro- 
scopic system B ("bath" or "reservoir"). Due to this 
contact the symmetries of the system S and the corre- 
sponding conservation laws are violated. If the contact 
between two systems is weak enough, the dynamics of 
interaction may be described in terms of transitions be- 
tween the eigenstates of a system S belonging to different 
irreducible representations of the group Gs generated by 



the operators which obey the algebra g s . If the operators 
describing transitions between these eigenstates together 
with generators of the group G5 form the enveloping al- 
gebra d 5 for the algebra g s , one may say that the sys- 
tem S possesses a dynamical symmetry characterized by 
some group D$. The dynamical symmetry group tech- 
nique offers mathematical tools for the unified approach 
to quantum objects, which allows one to consider not 
only the spectrum of a system iS, but also its response 
to external perturbation violating the symmetry G5 and 
various complex many-body effects characterizing the in- 
teraction between the system S and its environment B. 
In this paper we discuss a general algorithm of dynami- 
cal symmetry group approach to the few electron objects 
B + S and its practical application to the single-electron 
tunneling in complex quantum dots and single molecule 
transistors. This tunneling is described in a framework 
of the Anderson model with repulsive and attractive in- 
teraction between the confined electrons. 



II. HUBBARD OPERATORS GENERATING 
THE SPECTRUM OF NANOOBJECT 

Following the definition used in Ref. |Tl|, we define the 
dynamical symmetry group D 5 as a Lie group of finite di- 
mension characterized by the irreducible representations 
which act in the whole Hilbert space of eigenstates |^A) 
of the Schroedinger equation 

H\lX)=E l \lX) (2.1) 

describing the quantum system S. Here I is the index of 
irreducible representation and A enumerates the lines of 
this representation. The projection operators 

X$ = \l\){lii\ (2.2) 

play the central part in the procedure of construction of 
irreducible representations / of the group of Schroedinger 
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equation G s . The basic property of these operators is 
given by the equation 

X$\l'v) = 5 lv 8 tiV \lX) . (2.3) 

One may add to the set (|2.2[) the operators 

X^ = \lX)(l^\ (2.4) 

which project the states belonging to different irreducible 
representations (I ^ I') of the group G s one onto another. 
These operator may be also used for construction of the 
Lie algebras d 5 generating the spectrum of eigenstates of 
the Schroedinger equation and transitions between these 
states. Unifying the notations \IX) = A), we obtain the 
commutation relations 

[X AA ',H] = (E A ,-E A )H (2.5) 

The right hand side of Eq. (|2.5[) turns into zero provided 
the states A and A' belong to the same irreducible repre- 
sentation of the group G 5 . 

If a closed algebra d s exists for the set of operators 
(I2.2[) . (I2.4[) . then one may state that the system described 
by the Hamiltonian (12. ip possesses the dynamical sym- 
metry D s . This algebra is conditioned by the norm 

J2X AA = 1 (2.6) 

A 

and by the commutation relations for the operators 
AT AlA2 . In the general case these relations may be pre- 
sented in the following formal 

[X A ^,X A ° A % = I AiA M A2 a 3 t^ a ^ AiA4 (2.7) 

The "general case" implies that the Fock space includes 
the states which may belong to different charge sectors, 
i.e. changing the state Ai for the state A2 means chang- 
ing the number of fermions A/"a 2 — ^ A/ai in the many- 
particle system. If both M Al — Af\ 2 and Na 3 — A/a 4 are 
odd numbers (Fermi- type operators), the plus sign should 
be chosen in Eq. (I2.7[) . If at least one of these differences 
is zero or even number (Bose-type operators) , one should 
take the minus sign. 

The operators X AlA2 were exploited by J. Hubbard 
as a convenient tool for description of elementary excita- 
tions in strongly correlated electron systems (SCES). His 
seminal model of interacting electron motion in a narrow 
band, known now as the Hubbard modeler— was the 
first microscopic model of SCES for which the conven- 
tional perturbative approach based on the Landau Fermi 
liquid hypothesis turned out to fail. Now the realm of 
SCES is really vast, and the most of artificial nanostruc- 
tures in fact belong to this realm. In particular, complex 
quantum dots under strong Coulomb blockade are typical 
examples of short Hubbard chains or rings. 

The Hubbard operators (|2.4p obeying the commutation 
relation (|2.5I) provide a convenient tool for construction of 



the algebras generating the dynamical symmetry groups 
of the Schroedinger operator (H — E) or theresolvent op- 
erator R = (H — 2?) _1 . These operators may be used for 
construction of the generators of the appropriate group 
D s and irreducible tensor operators 0^ r > (scalars, r = 0, 
vectors, r = 1, and tensors r = 2,3...) which transform 
along the representations of the group D 5 : 

{ e r) =J2(MO ( g ) \A')X AA '. (2.8) 

AA' 

Here the index g stands for the components of irreducible 
tensor operator of the rank r. On the one hand, it is clear 
that the operators X AA generate all the eigenstates of the 
Hamiltonian H from any given initial state A'. The com- 
ponents of the operator form their own closed alge- 
bra, which characterizes the dynamical symmetry group 
provided the Hamiltonian H possesses such symmetry. 
Having in mind the application of this technique to the 
geometrically confined nanoobjects, we restrict ourself by 
the discrete eigenstates. 

The Clebsch-Gordan expansion (|2.8[) is the basic equa- 
tion which allows one to treat the dynamical symmetries 
of nanoobjects in a systematic way. The principal dif- 
ference between the dynamical symmetries of SCES and 
those of integrable models is that in the latter case the 
spectrum of the object and its dynamical symmetries are 
known exactly, while in the former case as a rule only 
some part of excitation spectra (usually its lower part) 
may be found analytically and classified by symmetry. 
This means that one may judge about the dynamical 
symmetry of the spectrum only within the definite en- 
ergy interval E. Respectively, the characteristic energy 
scale may be different for different problems. 

Our main subject in this paper is the Kondo effect in 
quantum dotsi 23 ' 24 The hierarchy of the energy scales in 
this problem is well know m 25 ' 26 The Kondo effect arises 
as a result of orthogonality catastrophe in the Anderson 
model,— where the conduction electrons in the Fermi sea 
of metallic electrodes play part of the subsystem B and 
the strongly correlated electrons in the quantum dot rep- 
resent the subsystem S. The largest energy scales in 
the Anderson model are the width of conduction band 
D in the subsystem B and the energy of Coulomb block- 
ade Q in the subsystem S. The electrons confined in 
the nanoobject (quantum dot) are characterized by the 
ionization energy ej. Next in the hierarchy of energies 
are the tunneling amplitude V and the tunneling rate 
r = irpoV 2 characterizing the process of electron tun- 
neling through the potential barrier, which separates two 
subsystems. Here po is the density of electron states at 
the Fermi level e F of the electron liquid in the leads. 
The Kondo effect arises in the single-electron tunneling 
regime under the restrictions of strong Coulomb block- 
ade Q. In this regime the charge transport between the 
source and drain electrodes constituting the subsystem B 
is realized as the electron cotunneling, where an electron 
from the source may tunnel into the dot S only provided 
another electron leaves the dot for the drain. Cotunnel- 
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ing, which arises in the fourth order in V , is characterized 
by the energy J. Finally, the energy scale of Kondo ef- 
fect Ek ~ s/DT cxp(— l/p J) characterizes the crossover 
from the weak coupling regime J -C 1 to the strong cou- 
pling, where J is enhanced due to the multiple creation 
of low-energy electron-hole pairs in the leads in the pro- 
cess of cotunneling. The Kondo energy also scales the 
excitations above the ground " Kondo-singlet" statei 25 ' 26 
The hierarchy of all these energies is 

L>,C/>e l >F>r>J>^ (2.9) 

An effective way to describe the crossover from the 
weak coupling to the strong coupling Kondo regime is 
the renormalization group (RG) approachi 2 ^— In this 
method the renormalization of parameters e,, T, J in (|2.9I) 
in the course of reduction of the energy scale £ from high 
energies ~ D,U to low energies still exceeding Ek is 
calculated. Our purpose is to describe this procedure 
in terms of dynamical symmetries which change in the 
course of reduction of the energy scale £. It was no- 
ticed that the multistage Kondo screening predetermines 
the non-universal features of the Kondo tunneling in the 
quantum dots with even occupation.— ~— In that case 
the relevant dynamical symmetry groups are SO{n) with 
n = 4 — 8.25i2& In this paper we will show that this lan- 
guage is useful already in the studies of the "ordinary" 
Kondo effect for quantum dots with odd electron occu- 
pation Af characterized by spin 1/2. The relevant Lie 
groups are SU (n) with n = 3,4. 



III. DYNAMICAL SYMMETRIES IN 
QUANTUM DOTS 

As was mentioned above, the dynamical symmetries 
of confined electrons in the quantum dot S are revealed 
in its interaction with the " Fermi bath" B of conduction 
electrons, The Anderson Hamiltonian describing the cou- 
pling between two subsystems reads 

H = H d + H b + H db (3.1) 

where three terms describe the nanoobject, the Fermi- 
bath and their coupling, respectively. The term Hdb in 
general case includes the direct coupling (quantum tun- 
neling of electrons between two subsystems), the direct 
interaction of Coulomb and exchange nature and the indi- 
rect (kinematic) interaction induced by the tunneling. If 
the symmetry of nanoobject is well defined, the Hamilto- 
nian H d may be diagonalized by means of projection op- 
erators (|2.2|) . and the generators of dynamical symmetry 
group (|2.8p arise in the interaction term Hdb in combina- 
tion with the operators describing the excitations in the 
Fermi bath. These symmetries cannot be treated in the 
same way as the symmetries of the integrable systems dis- 
cussed in the monographs MITT1 and the references therein, 
because the interaction not only activates the symmetry 



Ds of the nanoobject but also involves the charge, or- 
bital and spin degrees of freedom of the bath. This prin- 
cipal difference was pointed out in Refs. |H, where the 
quantum tunneling through an artificial molecule (dou- 
ble quantum dot) with even electron occupation M = 2 in 
presence of the many-particle interaction of Kondo type 
was described by means of the generators of the SO (4) 
group. 

To take the dynamical symmetries explicitly in the cal- 
culations of excitation spectra and in the studies of spin 
and charge transport in nanoobject, one should adhere 
the following paradigm: 37 

• When diagonalizing Hd use the projection opera- 
tors in accordance with Eqs. (|2.ip - (I2.5[k 

• Construct the operators X AA , which describe tran- 
sitions between all the states in the "supermulti- 
plet" of eigenstates of Hd belonging both to the 
same and to different irreducible representations of 
the symmetry group of the Hamiltonian H d and 
determine the relevant closed algebra generating 
the dynamical symmetry group Dg; 

• Rewrite H d b in terms of the configuration change 
operators (|2.4|) belonging to adjacent charge sectors 

• When projecting the original Anderson Hamilto- 
nian (13.11) on the subspace of low-energy states 
(A| . . . | A) by means of the Schrieffer- Wolff (SW) 
transformation 3 ^ or its generalizations, express the 
Hubbard operators which arise in this transforma- 
tion via the generators of corresponding dynamical 
symmetry group using expansion (|2.8p . 

To demonstrate this paradigm in action, let us con- 
sider the textbook example of a cell which may contain 
zero, one or two electrons with zero orbital moment. The 
Hamiltonian of this toy model 

H d = e d ^ d L d ia + Un ldt n ldi (3.2) 

is nothing but the single-site Hamiltonian describing the 
elementary cell of the non-degenerate Hubbard model 2 ^ 
with variable occupation number JV = 0, 1, 2 ("Hubbard 
atom"). Using definition (12.41) of the Hubbard operator, 
we rewrite H d in the universal form 

H d = Y,E A X AA (3.3) 

A 

where A = 0, a, 2 and the energy levels Ea are 

E =0, E f = E i = E 1 = e d , E 2 = 2e d + U. (3.4) 

It is convenient to arrange the energy levels in accor- 
dance with the available charge and spin sectors (Fig. 
QJi). The arrows connecting the levels Ea and Ea> corre- 
spond to the Hubbard operator A AA and its complex 
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FIG. 1: (a): Scheme of the energy levels for a Hubbard atom 
with the SU(4:) dynamical symmetry describing transitions 
between the states with occupation TV = 0, 1, 2. (b) The same 
for a reduced spectrum with the SU(3) dynamical symmetry 
describing transitions between the states with AT — or 2 and 
Af = 1. The Bose-like transitions with even 8Af = 0, ±2 are 
marked by the dashed arrows, the Fermi-like transitions with 
odd SAf = ±1 are marked by the solid arrows. 



conjugate. This scheme visualizes the Fermi-like and 
Bose-like operators (solid and dashed lines, respectively), 
which obey the commutation relations (|2.7I) . There are 15 
such operators forming a closed superalgebra containing 
both commutators and anticommutators. The Hubbard 
operators may be regrouped into the generators of the 
SU(4) group, known as the Gell-Mann matrices of 4th 
rank Ai — A15 (see Appendix). Thus the generic dynam- 
ical symmetry of Hubbard atom which is realized within 
the energy interval £ ~ U, D is SU(4). 

Reduction of the energy scale to the interval < £ <C 
U results in quenching the doubly occupied levels and 
corresponding reduction of the dynamical symmetry from 
SU(4) to SU(3) (Fig. [TJd). The algebra generating this 
group contains eight Gell-Mann matrices of the 3rd rank 
Ai — As and the same number of Hubbard operators. Re- 
lations between the matrix of Hubbard operators and 
the Gell-Mann matrices for this group are also presented 
in Appendix. Further reduction of the energy interval 
£ <C £j results in complete suppression of charged sectors 
Af 7^ 1, so that we are left only with spin states a =t,l- 
In this limit the dynamical symmetry is the same as the 
symmetry of the Hubbard atom, and the corresponding 
Lie group is SU(2). 

Mathematically, non-trivial dynamical symmetries are 
described by semisimple groups, possessing ideals. These 
ideals are formed by some triads of Gell-Mann matrices, 
e.g. matrices Ai, A2, A3, which form a subgroup SU{2) of 
a group SU(n). If the states in the Fock space for the 



Hubbard atom are ordered as 

$4=(t I 2) (3.5) 

then the first three Gell-Mann matrices \\ — A3 are re- 
lated to the spin states in the charge sector Af = 1 (see 
Appendix) . 

It is expedient to rewrite the original Hamiltonian (|3.3[) 
in terms of the generators of the group SU(A) in the case 
where all four eigenstates (|3.4p shown in FigQJi are taken 
into account, and in terms of the SU(3) generators in 
the case when the polar states with Af — 2 are frozen out 
(Fig [lb). In the full space $4 we obtain by means of (|A3|) 

rrSU(4) _ E / 4 \ 

Ei ( 2 2 \ 

+ jl.+ ff 1 -v& 15 ) (3-6) 

Here the notation X p is used for the Gell-Mann matrices 
X p defined in the Fock space (|3.5[) . The Zeeman term hS z 
acting in the charge sector Af = 1 is also added. In the 
reduced Fock subspace 

*3 = (t I 0) or (t I 2) (3.7) 

the Hamiltonian of the Hubbard atom rewritten with the 
use of Eqs. (IA4[) acquires quite compact form 

h! U(3) = f (l - V3X 8 ) + I (l + V3A 8 ) + \X,. 

(3.8) 

The Hubbard atom is a minimal model which can be 
used for description of a quantum dot with variable occu- 
pation Af coupled with the bath by means of the tunnel- 
ing channel. The equilibrium occupation of the dot may 
be changed by means of injection of an electron or a hole 
from the metallic reservoir J£ This occupation fluctuates 
dynamically due to the single electron tunneling (SET) 
between the dot and the leads. The Coulomb blockade 
parameter Q plays the same part as the Coulomb repul- 
sion U in the original Hubbard model. In the general 
case of, say, planar quantum dot the energy spectrum 
of a quantum dot contains many discrete states without 
definite angular symmetry. Only the highest occupied 
(HO) and the lowest unoccupied (LU) states are involved 
in single electron tunneling through such quantum dot. 
The Hamiltonian of subsystem S in the Hamiltonian (|3.1[) 
has the form 

Hd = J2 e i4<r d i° + ^int + Q ("dot " (3.9) 

3 

Here the index j enumerates the levels bottom-up. Hint 
is the electron-electron interaction in the quantum dot, 
Usually the self-consistent Hartrcc term is included in the 
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definition of discrete levels Sj , and the relevant contribu- 
tion to Hi n t is the exchange between the electrons occupy- 
ing different levels of a neutral quantum dot. Q = e 2 /2C 

is the capacitive energy of the dot, ndot = Ylja^ ^ja^ja- ^ s 
the number of extra electrons or holes which are injected 
in the dot due to tunneling described by the Hamiltonian 
H db . 



l—s,d jk<r 



H.c. 



(3.10) 



Corrections to the capacitive energy take into account the 
capacitance of the gate C g and the gate voltage v g . If the 
hierarchy of the energy scales 

Q>(fc,J)>Wy (3.11) 

takes place (Se is the interlevel spacing between the HO 
and LU states, J is the exchange coupling constant), then 
one may assert that the charge transfer through the quan- 
tum dot occurs in the SET regime. 






FIG. 2: Upper panel: Variation of the energy of the quantum 
dot E e \(N) as a function of the gate voltage. Lower panel: 
corresponding variation of addition energies for electron and 
hole excitations relative to the Fermi level in the leads. See 
the text for further discussion. 

Variation of the energy spectrum and the occupation 
of the quantum dot as a function of a gate voltage is 
exemplified in Fig. [3J The "Hubbard parabolas"^— 
represent the energy E e \(M) of the isolated quantum dot 
with the Hamiltonian (13.9[) . Three subsequent diagrams 
for the occupation Af — 1 show the asymmetric configu- 
rations with quenched zero and two electron occupation 
(the side diagrams) and the configuration with particle- 
hole symmetry (the middle diagram). The single parti- 
cle excitations are the addition and extraction energies 
E(Af) — E(Af =F 1) which should be compared with the 
chemical potential of the bath (the Fermi energy) - see 
the lower panel of Fig. [2] 

Following the above scheme, we express the Hubbard 
Hamiltonian (|3.6[) and the tunneling term Hdb 



H, 



(ib 



(t k d 



1 r 

ka^ka 



H.c. 



(3.12) 



via the operators from the triads (| A2|) which are con- 
nected with the original Hubbard operators X AA acting 
in the space (|3.5p by the following relations: 

T+ = X n , J- = X l \ T z = X n - X u 

v+ = v t0 , v-=x ot , v z = x n -x 00 
u+ = x^ u- = x°\ \J Z = X U -X 00 
w+ = x n , v\r = x 2t , \n z = x n - x 22 

Y + = x i2 , Y~ = X 2i , Y z = X u - X 22 

Z+ = X 02 , Z-=X 20 , Z Z =X 00 -X 22 (3.13) 



Equations (|3.13l) realize the general expansion scheme 
(|2.8p for the irreducible vector operators in the group 
SU(4:). The triad T is nothing but the set of spin 1/2 op- 
erators (S + , S~ , 2S Zl ) acting in the charge sector N = 1. 
The triad Z describes the two-particle excitations (N — 
-H> Af = 2). The rest four triads describe transitions 
between different charge sectors (Af = 1 <-> J\f = 0, 2). 
These operators enter the Anderson Hamiltonian corre- 
sponding to the Hubbard parabolas of Fig. [5J 

The Hamiltonian JJ^ U may be expressed via the z- 
components of irreducible vectors (|3.13l) by means of the 
following relations 



P = V z + U z = X 11 - 2X Q0 
Q = W z + Y z = V 11 - 2X 2 ' 
R = P - Q = 2(X 22 - X 00 ) 



(3.14) 



and the completeness condition (|2.6I) which in this case 
reads 

x 00 + x 11 + x 22 = i, x 11 = J2 xaa ( 3 - 15 ) 



[see also Eq. (|A"5l) ]. We find from Eqs. (I3TT4]) and (j3"T5j) : 
X°° = i-i(3P-Q) 
V 22 = i + i(P-3Q) 



X" = i + i(P + Q). 



(3.16) 



Then the general SU (4) configurations (the first and the 
third parabolas in Fig. [2]) are described by the Hamilto- 
nian 

jjsu(4) = 2E i +E + E 2 1+ h j 

Ell.p + Ell.Q + ^.R (3.17) 

4 4 8 

where 1 is the unit matrix in the Fock space $4, By = 
Ei — Ej are the addition/extraction energies. Thus the 
operators P/4, Q/4, R/8 and T z /2 describe all Fermi- and 
Bose-like excitations shown in Fig. QJi. In the degenerate 
case Eq — E2 = E e (second parabola in Fig. [2]) this 
Hamiltonian reduces to 

H s d ^^^±^-l + h --T z + ^-(Q + P) (3.18) 
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The tunneling term H^b also may be expressed via the 
generators of SU(4) group, namely via the ladder opera- 
tors: 

H S db U[i) =^t fc (V t +Wt) Cfet + (Ut-Y-) Cfc ,+H.c. (3.19) 



In the strongly asymmetric situations (the side configu- 
rations in Fig. [3J where the excitation Egi is soft, whereas 
the excitation E21 is frozen out or v. v.), the symmetry of 
the dot is reduced from SU(4) to SU(3). Respectively, 
the system (|3.16[) reduces to 



1 P 



2 P 



X -3~3> X ~3 + 3' 
or, in terms of operators U, V 



(3.20) 



^ 1 2V 2 - u 2 , , 1 2U 2 - v 2 , 

x n = _ + — ! i x u = - + — -.(3.21) 

3 3 3 3 v ; 

The Anderson Hamiltonian acting in the space $3 has 
the form 

Hf (3) = 2 -^±3> . 1 + |£ . ( Uz + V z ) + ^. T 2 (3.22) 



H. 



SU{3) 
db 



*k [( V+C fct + U+c H ) + He] . (3.23) 



Thus the operators describing the charge Hubbard exci- 
tations in SU(3) subspace $3 (|3.T[) are U, V, whereas the 
spin excitations are described by the conventional spin 
operator S = T/2. 

The dynamics of charge and spin excitations in this 
case is predetermined by the commutation relations for 
the group generators. The operators belonging to the 
same subgroup (triad) commute in accordance with the 
standard SU(2) relations 



[0 2 ,0 ± ]=±20 ± , [0+,0-] = 2 



(3.24) 



The non-zero commutation relations between the opera- 
tors belonging to different triads ensure complex dynam- 
ical properties of Hubbard-like SCES. 

[U±V*] =±T=f, [U±,V Z ] = T U ± , 
[U z y ± }=±V ± , [U 2 ,V 2 ] =0. (3.25) 

Respectively, the non-zero anticommutation relations are 

2 + V 2 -2U Z 



{u+,u-} 
{v+,v-} 



2 + U 2 - 2V Z 



(3.26) 



Then the excitations in the charge sector are described 
by the Green functions, which may be found directly from 
equations of motion for the generators of SU(3) group. 

G v = «V-(t)V+(0))), G„ = ((ir(t)U + (0))), (3.27) 



Respectively, the excitations in the spin sector are given 
by the Green functions 



G a = «S-(f)-S+(0)>, 



(3.28) 



Here the double brackets stand for thermal averaging and 
time-ordering operations specified for the retarded, ad- 
vanced or causal Green function. These functions can 
be easily found in the atomic limit where only the term 
H^ U is retained. Solving equation of motion for the 
"Fermi-like" Green functions which describe excitations 
in the charge sector, one gets by means of the commuta- 
tion and anticommutation relations (I3.25|) and (|3 . 26|) : 



G v (u) 



% (2 + (V 2 ) -2(U 2 ))/3 
2n uj — ed 

1 (2+(U 2 )-2(V 2 ))/3 



2vr 



(3.29) 



Using the definitions (|3. 201) and (|3.2ip . we see that the 
numerators in the Green functions (l3~29)) are nothing but 
the averages (X 00 ) + (X n ) and (X 00 ) + (X u ), so that 
these functions are indeed the atomic Green functions for 
the Hubbard model 22 rewritten in terms of the generators 
of the SU(3) group. 

The " Bose-like" Green function G s which describes the 
excitations in the spin sector in theatomic limit has the 
usual form 



G. = 



2ir uj 



(3.30) 



A. Three-fold way for the Hubbard atom 

As was noticed in the sixties, - various families of 
hadrons are classified in accordance with the irreducible 
representations of SU(3) group (see also Ref. HI). In 
particular, 18 baryons form two multiplets corresponding 
to representations D^ 11 ^ (the octet of baryons with spin 
1/2) and Z)( 30 ) (the decuplet of baryons with spin 3/2). 
The octet of spinless mesons also transforms along the 
representation D^ n \ The higher representations of the 
SU(3) group are realized in the physics of strong interac- 
tion because these " composite" particles possess not only 
the spin and charge but also the isospin and hypercharge 
quantum numbers, and the SU(3) symmetry character- 
izes the latter variables. The elementary particles obey- 
ing the SU(3) symmetry are the colored quarks. The 
SU (3) symmetry in the hadron multiplets under strong 
interaction is satisfied only approximately due to exis- 
tence of electro-weak interaction, so that this symmetry 
may be treated as a dynamical symmetry in the original 
sense of this notion. 

The Hubbard atom with frozen doubly occupied states 
possesses only two quantum numbers, namely spin and 
charge. Therefore, the multiplet of Hubbard states is de- 
scribed by the lowest irreducible representation Z?( 10 ' of 
the SU (3) group. To construct this representation, one 
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should recollect that the two of eight Gell-Mann matrices 
can be diagonalized simultaneously. Following Ref. H^, 
we choose the representation with diagonal matrices T z 
and Q. Then the set of allowed states is defined by two in- 
teger numbers A, \i so that the eigenstates are determined 
as 

M T = A + M , Mq = ^(A — fx). (3.31) 

The whole set of eigenstates form a two-dimensional tri- 
angular lattice on the plane (M T , Mq). Each irreducible 

representation D x ^ is marked by the indices A, ji corre- 
sponding to the state with the maximum eigenvalue Mq 
and the maximum value of M T possible at this Mq . Then 
the rest states forming this irreducible representation are 
constructed by means of the ladder operators T^, l^, 
acting on the state \Mq, M T ). 

This procedure results in construction of the stars of 
basis vectors D x ^ and the polygons connecting the points 
generated by the ladder operators subsequently acting 
on the point (Mq, M T ). In the case of baryon family 
the corresponding multiplets are the hexagon with dou- 
bly degenerate cental point for representation Z?( n ) and 
the triangle with ten point in its vertices and on its sides 
for representation D^ 30 \ In the case of Hubbard atom 
the multiplet is represented by a triangle (Fig. la- 
beled in accordance with the state with the highest quan- 
tum numbers A = 1, fx = 0, which corresponds to the 
state |jV, cr) = |l,f) of the Hubbard atom. Two remain- 



are presented in the following table 



(-1,1/3) 



-1 V 



T- d,l/3) 



/I 

v- 

(0,-2/3) 



M 7 



FIG. 3: Irreducible representation D' 10 ' for the set $3 

ing components of the multiplet <I>3 may be generated 
from the state |l,t) by means of the ladder operators 
T = X^ and V~ = X ^. First of these operators corre- 
sponds to the "Bose-like" excitation with spin 1, and the 
second one is the "Fermi-like" excitation with spin 1/2. 
The triangle D^ 10 - 1 is closed by means of the operator 
U + = X^°. The interrelations between the values of the 
parameters A, [i, the eigenvalues of the operators T z and 
Q, and the eigenvalues |A) of the Hubbard Hamiltonian 



A 


/'■ 


M T 


Mq 


A 


1 





1 


1/3 


u 





-1 


-1 


1/3 


d 


-1 


1 





-2/3 


h 



(3.32) 



Here the notations u, d, h are used for the spin up, spin 
down and hole states, respectively. 

Thus we see that the dual nature of the Hubbard oper- 
ators manifested in the superalgebra with the commu- 
tation relations (|2.7I) allows one to use them for con- 
struction of the generic su(3) algebra formed by the spin 
and pseudospin operators with the commutation relations 

Like in the case of baryons and mesons, this symmetry 
is violated due to interaction with other subsystems. In 
our case this is the Fermi bath B. The source of this inter- 
action is the tunneling coupling given by the Hamiltonian 
fifth ■ In analogy with the term offered by Gell-Mann 
and Ne'eman for the multiplet of light hadrons one may 
use the term "three- fold" (ternary) way for the SCES 
with approximate SU(3) symmetry. Generalization of 
this description for the SU(4) group is straightforward. 
In this case the phase space for the irreducible represen- 
tations is defined by the eigenvalues of the operators P, 
Q, T z , and the lowest irreducible representation of this 
group Z)( 100 ) is represented by a triangular pyramid in 
this 3D space. 

In the next section we will see how the hierarchy of dy- 
namical symmetries of the Hubbard atom manifests itself 
in the RG evolution of the Anderson-Kondo problem. 



IV. TWO-STAGE RENORMALIZATION GROUP 

FOR 517(3) AND 517(4) ANDERSON 
HAMILTONIAN 

The RG method is based on the idea of renormaliza- 
tion of model parameters, which are relevant at low en- 
ergy as a result of the change of the scale of high energy 
excitations^ If the model is renormalizable, any such 
parameter P(e) may be represented as 



P(e)-P[(l + «)e] 



neP'(e) 



(4.1) 



where k is positive infinitesimal and the prime stands for 
the derivative. The quantity —KsP'(e) is the contribu- 
tion to P(s) from the high-energy states which are to be 
integrated out, preserving the form of P(e) but changing 
its scale. Adopting this approach, we immediately notice 
the inevitability of the three- or two-stage RG procedure 
as a direct consequence of several energy scales inherent 
in the Anderson model and the dynamical SU(n) sym- 
metry of its excitation spectrum with n — 4 or 3. 

Taking as an example the first of three Hubbard 
parabolas in Fig. [21 we see that in this case the high- 
est energy scale is the addition energy £ ~ £rf + Q — sp. 
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The corresponding generators of SU (4) group are W, Y. 
The next energy scale E ~ e F — e d is the extraction en- 
ergy, and the relevant generators are U, V. The lowest 
energy scale £ <C t 2 /ed is introduced by the second-order 
cotunneling processes from the dot to the leads, which are 
accompanied by the spin flips in the dot and creation of 
the low-energy electron- hole pairs in the leads. The vec- 
tor operator T is responsible for these processes. In other 
words, we arrive to the the renowned Jefferson-Haldane- 
Anderson renormalization group (RG) procedure . 31 ' 32 
Basing on the symmetry analysis of preceding section, 
the RG procedure may be described in terms of the gen- 
erators of SU (4) group and its subgroups with reduction 
of the symmetry 517(4) SU(3) -> SU{2) following the 
reduction of the energy scale £. 

Let us rederive for example the scaling equations for 
the two stage RG SU(3) — > SU(2) realized in the limit 
U — > oo in these new terms. In the Anderson model 
the change of the energy scale in Eq. (|4.ip means the 
contraction of the electron bandwidth D — > D — 5D in Hjj, 
The renormalized quantities are the self energies of 
the Green functions (|3~2"7f and (pHSl) . T) = u,v,s. The 
tunneling Hamiltonian (|3.23[) gives the second-order self- 
energy part for the Green functions G v , G u (|3.29l) 



e<j — Ew H — 



de 



E 



10 



(4.2) 



where 



•E 

k 



\t k \ 2 8{s F -S k ) 



is the spin-independent tunneling rate. The transfor- 
mation (|4.ip results in the Jefferson-Haldane scaling 
equatio: 



,31,32 



de d _ _T_ 
dD ~ ttD 

with the scaling invariant 



it I r 



(4.3) 



(4.4) 



Thus, the evolution of the resonance level is determined 
by the vectors U, V operating in the charge subsectors of 
the group SU(3). The same second-order processes gen- 



erate the four-tail vertices ~ V+U cLc fc ,,p U + V c 



"k^k' 



etc. Using the commutation relations (I3.25[) . these ver- 
tices are combined in the conventional Schrieffer- Wolff 
exchange interaction 



(4.5) 



where s = N^ 1 ^2 kk , c^, (T fc fc , cr , is the local spin opera- 
tor for conduction electrons, f is the the vector of Pauli 
matrices J| ~ t kp \ 2 / E w is the indirect Kondo exchange. 
The scaling equations for this Hamiltonian may be de- 
rived by means of the Anderson's RG procedure. 



In the symmetric configuration (middle parabola in 
Fig. [2J the Jefferson - Haldane - Anderson scaling theory 
is in fact the manifestation of reduction of the dynamical 
symmetry SU{A) — > SU(2), in which the charge excita- 
tions represented by the vectors (W,Y, U,V) are frozen 
out in the process of renormalization and the subgroup T 
describing the spin degrees of freedom in the charge sec- 
tor Af = 1 represents the low-energy part of the spectrum 
responsible for the Kondo singularities. Five vectors of 
six triads available in the Gell-Mann set are involved in 
this two-stage procedure. 

Another choice of 15 linearly independent generators 
of the ST/ (4) dynamical symmetry is possible in the case 
where instead of the spin degeneracy of the ground state 
with Af = 1 the charge degeneracy of the two singlets 
with Af — 0, 2 is realized. Such a possibility arises in 
the negative U Anderson model42r— In this configura- 
tion the vector T is excluded from the renormalization 
procedure due to quenching of the sector Af = 1 at 
low energies. Instead the vector Z is involved in for- 
mation of the Kondo singularities. The attractive in- 
teraction between the electrons in the nanoobject in this 
model stems from the strong electron-phonon interaction. 
Starting with the Anderson - Holstein model where the 
phonon subsystem is represented by the single Einstein 
mode with the energy 51o, one may perform the canonical 
transformation, 43 which transforms the electron-phonon 
interaction into the polaron dressing exponent for the 
electron tunneling rate, the polaron shift of discrete elec- 
tron levels and the phonon mediated electron-electron in- 
teraction. The latter renormalizes the Hubbard interac- 
tion term in the Anderson Hamiltonian 



u' = u- 2\ 2 n a 



(4.6) 



Here A is the electron-phonon coupling constant. In the 
limit of strong electron-phonon coupling the energy gain 
due to the phonon mediated interaction overcomes the 
energy loss due to the Hubbard repulsion, and one comes 
to the case U' < 0. The negative U model may be real- 
ized in the single electron molecular transistors The 
interaction (14.61) should be included in the term Hd, so 
that in the negative U case the Hubbard parabolas for 
the energy spectrum are reversed relative to the usual 
shape shown in the middle configuration of Fig. [5] The 
"turned over" diagrams corresponding to the two nearly 
symmetric configurations shown in Fig. |4j Like in the 
positive U case, the transitions between the levels in the 
Hubbard supermultiplet are described by the operators 
(|3.13[) generating the Sf/(4) dynamical symmetry group. 
We consider here the configurations, where the singlet 
states | A) = |0), |2) are degenerate or nearly degenerate, 
and the spin doublet |A) = | f), | 1) is an excited virtual 
state in the cotunneling processes. The two configura- 
tions presented in Fig. [4] correspond to the empty and 
completely filled two-electron shell of the Hubbard atom. 
They are connected by the particle-hole symmetry trans- 
formation, so it is enough to discuss one of them. 
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E, n /2 



FIG. 4: Upper panel: inverted Hubbard parabolas for the 
negative U Hubbard atom in the cases of empty and doubly 
occupied shells. The interlevel transitions are described by 
the operators generating the S'C/(4) dynamical group. Lower 
panel: single-electron levels corresponding to the transitions 
shown by the arrows in the upper panel (see the text for fur- 
ther explanation). 



We will show below that the negative U Anderson 
model may be formally mapped on the positive U model, 
by means of the multistage RG method, which general- 
izes the Jefferson - Haldane - Anderson procedur o 28 i 31 ' 32 
mentioned above. In the positive U case after freez- 
ing out the high-energy excitations E 01 and E21 corre- 
sponding to injection of a hole or of an electron into the 
singly occupied quantum dot at the Jefferson-Haldane 
stage of the renormalization, one arrives at the Ander- 
son stage of Kondo screening of spin excitations in the 
sector J\f = 1 described by the vector operator T. In the 
negative U model the spin excitations are exponentially 
suppressed from the very beginning. After freezing out 
the charge excitations E w and E 12 generated by the op- 
erators U,V,W,Y we are left only with the two-particle 
charge excitations E 20 generated by the operator Z. 

Since the Jefferson-Haldane stage of the RG procedure 
is realized exactly in the same way as in the powitive 
U Anderson model, we concentrate on the second stage, 
where the SU(4) dynamical symmetry group is reduced 
to its SU(2) subgroup represented by the triad Z. These 
operators act in the subspace $2 = (0,2). The effective 
SW Hamiltonian in this subspace reads 



Jx 



tfcotun = (Z+B- 



Z~B H 



NJ\\Z Z B Z 



(4.7) 



where the components of the vector Z are presented in 
the last line of the system (|3.13|) ]. The components of 
the vector B defined in the space of two-particle itinerant 



excitations are 

B+ = N-^iA'V B-=N- 1 J2 c uCk'V 



B r = N~ 



(4.8) 



kk' <t 



These operators obey the su(2) commutation relations 

[B+B-] = B„ [B Z ,B ± ]=±B ± (4.9) 

The transversal part of the Hamiltonian (|4.7[) describes 
the tunneling of singlet electron pairs between the leads 
and the molecule, whereas its longitudinal part stems 
from the band electron scattering on the charge fluctua- 
tions. 

Thus the Hamiltonian of two-electron tunneling 
is formally mapped onto the anisotropic Kondo 
HamiltonianiLi 2 .^ (see FigEJ). The origin of this 
anisotropy is the polaron dressing of tunneling matrix 
elements^ This dressing is different for the two-electron 
cotunneling and the electron scattering coupling param- 
eters in the Hamiltonian (|4.7I) . In the strong electron- 
phonon coupling limit, (A/f^o) 2 = S ^> 1 



^ = (2|0) 



3 -2(A/n ) 



J\ 



(4.10) 



The eventual source of this anisotropy is the overlap be- 
tween the phonon wave functions for a molecule in the 
charge states N = and Af = 2, i.e. the Huang-Rhys 
factor S. In a framework of the Anderson RG scalng pro- 
cedure this means that the renormalization diagrams for 
the two models are the same, namely the diagrams in the 
first and the second columns of Fig. [5] are mapped on the 
longitudinal and transversal components of the Kondo 
exchange Hamiltonian. (|4.5[) The mapping procedure im- 
plies the substitution S —> Z, s — > B. The scaling equa- 
tions, which follow from these equations are the same as 
for the conventional anisotropic Kondo model,^ namely 



dj\\ 
drj 



•9 

-3±, 



drj 



■J±J\\ 



(4.11) 



(ji = poJi). In the case of strong anisotropy (14. 10)) solu- 
tion of this system gives for the Kondo temperature the 
following equatio^ 1 



n 



Z?exp 



"2F [jh 



(4.12) 



The last equation in (I4.12[) is valid in the limit of strong 
electron-phonon coupling S > 1. Generally, the polaron 
narrowing of the tunneling rate results in a noticeable 
decrease of Tx in comparison with its value for the con- 
ventional Kondo effect. 
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FIG. 5: RG diagrams in the space $2 = (0,2). Upper panel: 
the bare vertices J± for the two-electron tunneling and Jy 
for the charge scattering. Lower panel: the diagrams for 
the second-order renormalization of these vertices. Solid lines 
stand for the conduction electron states, dashed lines denote 
the charge states of the molecule. 



V. CONCLUDING REMARKS 

In this paper we have shown that the "atomic" part Hd 
of the Anderson Hamiltonian may be treated as a four- 
level system, so that the interlevel transitions induced 
by the tunneling term H c ib activate the implicit SU(4) 
dynamical symmetry of the model. The symmetry is the 
same for the negative and positive U models. However, in 
spite of the formal similarity between the effective Hamil- 
tonians for the single electron cotunneling and the elec- 
tron pair cotunneling, the background physics is different 
in two versions of the Anderson model. In the positive 
U Anderson model the tunneling in the middle of the 
Coulomb window arises exclusively due to the many-body 
Abrikosov - Suhl resonance. In the negative U model 
the resonance conditions for the two-electron tunneling 



arise at E02 = irrelative to the many body particle-hole 
screening mechanism, so that the zero bias anomaly in the 
tunneling conductance exists already at T ^ Tk, as well 
as the finite bias anomaly at £^02 7^ 042r— One may say 
that the Anderson orthogonality catastrophe^ respon- 
sible for the many-body Kondo-like screening at low T 
only enhances the two-electron tunneling resonance al- 
ready sharpened due to non-orthogonality of the phonon 
clouds measured by the Huang - Rhys factor (|4.10|) . The 
finite difference £02 7^ in the negative U model is equiv- 
alent to the finite magnetic field in the positive U model: 
it results in the appearance of two split finite bias peaks 
in the tunneling conductance. 

Having in mind all these differences, one may state that 
the multistage RG procedure reveals the hierarchy of re- 
duced dynamical symmetries SU(A) -> SU(3) -> SU(2) 
in the Anderson model both with the Hubbard repul- 
sion for odd occupation and with the Hubbard attraction 
for even occupation. In this paper we confined ourselves 
with rephrazing the results already known in the Kondo 
physics, but we believe that the use of operators generat- 
ing the set of the eigenstates of the zero Hamiltonian in 
the Anderson and Hubbard models and thus revealing its 
internal symmetry may facilitate the description of vari- 
ous dynamical properties of strongly correlated electron 
systems. In particular, new bosonization and fcrmioniza- 
tion procedures for the generators (|3.13p alternative to 
the standard representations for the SU(n) groups^— 
may be elaborated. We will turn to these procedures in 
the forthcoming publications. 



Appendix A: GELL-MANN MATRICES AND 
HUBBARD OPERATORS 

Here we summarize for the sake of convenience some 
properties of the Gcll-Mann matrices of 4th rank and 
their realization in the Hubbard and Anderson models. 
The canonical form of these matrices describing the sym- 
metry of four-level systems is 
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10 
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1 

71 



10 
10 
V -3 / 



(Al) 



(see, e.g^S). First eight matrices contain the 3-rd rank 
Gell-Mann operators of the SU (3) group as submatrices. 
The operators Ag — A15 generate transitions between the 
triplet and the fourth level. 

One may construct a subgroup SU{2) of the group 
SU{n) for any 2D subspace of the effective Fock space. 
There are three such "triads" grouped in three vectors 
T, U,V for the group SU(3) with the symmetry opera- 
tions acting in the 3D space $3 (cf. Ref. l39h - Adding 
fourth dimension provides three more vectors W, Y, Z rep- 
resenting the generators of the group SU (4) together with 
the first three vectors. 

T ± = i(A 1 ±iA 2 ), J Z = X 3 

U± = i(A 6 ± iX 7 ), U 2 = i(-Aa + VSX S ) 

V± = i(A 4 ± i\ 5 ), V z = i(Aa + V3X S ), 



= i (A 9 ± iX w ) , W 2 - X - f A3 + -J=A 8 + Aa 18 
i (An ± iAia) , Y z = 1 (-A 3 + -^=A 8 + -i=Ai 5 



- (Ai 3 ± iX u ) , Z z = — (-A 8 + V2X 



(A2) 



The original Hubbard operators A AA are represented 
via the Gell-Mann matrices in the basis $4 (|3.5[) in the 



following way: 

Xto = i(A 4 + iA 5 ), A°t = i(A 4 -*A 5 ), 
X^ = i(A 6 + iA 7 ), A°+ = i(A 6 -*A 7 )/2, 
^ 2t = i(A 9 -*Ai ), At 2 = I(A 9 + zAi ) 
A 2 + = i(A u -iA 12 ), A^ 2 = i(An+zAi 2 ) 
X n = I(A 1 +iA 2 ), X^ = i(A 1 -iA 2 ) 



A 
A 00 



^(Ai3-»Ai 4 ), A 02 = 7j(Ai3 + iXu) 



1 + 2A 3 
1-2A 3 



2 

7^ 
2 

7^ 



2 
2 

76 



A15 

Al 5 



u 1_ 7! A8 + ^ Al5 



X 22 = i(l-V6A«) 



(A3) 



The tree operators T 1 * 1 , T z from the first triad in the set 
(|A2I) describe the spin-flip excitations in the homopolar 
subspace J\f — 1 of the Hubbard atom. The three opera- 
tors Z , Z 2 from the last triad may be used in the descrip- 
tion of excitations in the two-particle sector J\f = {0, 2} of 
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the Hubbard and Anderson model. The operators form- 
ing the triads U and V intermix the states from the charge 
sectors Af — and Af = 1, and the operators W and Y 
do the same for the sectors N = 2, W = 1. 

In many physical applications the reduced Anderson 
and Hubbard Hamiltonians with U — > oo are exploited. 
In this limit the doubly occupied state |2) is completely 
suppressed. In the appropriately reduced Fock space $3 
(I3.7[) possessing the SU(3) symmetry the system (|A3I) 
transforms into 



x n 
x n 

A 00 



A ot = 




- «A 5 ), 


x oi = 


i(A« 


-iA r )/2, 


X** = 




-iA 2 ) 













1 f 2 A 1 A 



i C 1 - \/3A f 



Within each triad the standard Pauli commutation re- 
lations (|3.24[) for the components are valid. The com- 
mutation relations between the operators from different 
subgroups are described by more complicated structure 
factors^, These relations in our case may be derived from 
the general commutation relations (|2.7p for the Hubbard 
operators (see the main text). 

Two matrices used in the irreducible representation of 
the SU(3) group are 




(A5) 



/ 1 

T 2 = 0-10 
y 

Their eigenvalues are 

M T = 1,-1,0; M Q = 1/3, 1/3,-2/3 (A6) 

for the states | f), | 1), |0), respectively. 
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